Inversion of Sequence of Diffusion and Density Anomalies in Core-Softened Systems 



O 

(N 

oo 



o 

C3 



I 

c 

O 

o 



> 
o 



X 



Yu. D. Fomin, E. N. Tsiok, and V. N. Ryzhov 
Institute for High Pressure Physics, Russian Academy of Sciences, Troitsk 142190, Moscow Region, Russia 

(Dated: September 9, 2011) 

In this paper we present a simulation study of water-like anomalies in core-softened system in- 
troduced in our previous publications. We investigate the anomalous regions for a system with the 
same functional form of the potential but with different parameters and show that the order of 
the region of anomalous diffusion and the region of density anomaly is inverted with increasing the 
width of the repulsive shoulder. 
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I. INTRODUCTION 

It is well known that some liquids demonstrate anoma- 
lous behavior in some regions of thermodynamic pa- 
rameters. The most common and well known exam- 
ple is water. The water phase diagrams have regions 
where a thermal expansion coefficient is negative (den- 
sity anomaly), self-diffusivity increases upon compres- 
sion (diffusion anomaly), and the structural order of 
the system decreases with increasing pressure (struc- 
tural anomaly) Later on it was discovered that 
many other substances also demonstrate similar behav- 
ior. Some typical examples are silica, silicon, phosphorus 
and many others |5j-|20j. 

It is reasonable to relate this kind of behavior to the 
orientational anisotropy of the potentials, however, a 
number of studies demonstrate water-like anomalies in 
fluids that interact through spherically symmetric core- 
softening potentials with two length scales (see, for ex- 
ample, the reviews |2l| and [22j|). A lot of different core- 
softened potentials were introduced. All such systems 
can be approximately divided into two classes: shoulder- 
like, composed of a hard-core and a repulsive shoulder, 
softening the core, and ramp-like systems, composed of a 
hard-core and a repulsive ramp, establishing two compet- 



ing equilibrium distances [2ll . |22J . Anomalous behavior 
usually takes place in ramp-like systems, while it is not 
observed in purely shoulder-like ones. However, if in- 
stead of "hard" step one considers a smoothed one the 
anomalies appear (see, for example, |2ll424 |57 H61 |. 

As it was found in experiments [l| and simulations 
[H, IH , the water anomalies have a well-definite sequence: 
the regions where these anomalies take place form nested 
domains in the density-temperature [3| (or pressure- 
temperature Q) planes: the density anomaly region 
is inside the diffusion anomaly domain, and both of 
these anomalous regions are inside a broader structurally 
anomalous region. This water-like behavior was found 
in systems with core-softening potentials (2lT[23l [25l - l6l1 | . 
However, in other anomalous systems the sequence of 
anomalies may be different. For example, the hierarchy 
of anomalies for silica is different compared to water. In 
this case the diffusion anomaly region contains the struc- 



tural anomalous region which, in turn, has the density 
anomaly region inside [H^]. To our knowledge this is the 
only example of such inversion of the order of the anoma- 
lies. 

This paper presents a simulation study of anomalies in 
core-softened system introduced in our previous publica- 
tions 57-6l|. We investigate the anomalous regions for 
a system with the same functional form of the potential 
but with different parameters and show that the order 
of the region of anomalous diffusion and the region of 
density anomaly is inverted. 

The article is organized as following: Section II 
presents the system and methods, Section III describes 
the results and gives theirs discussion and Section IV 
contains conclusions. 



II. SYSTEMS AND METHODS 

The system we study in the present simulations is 
Smooth Repulsive Shoulder System (SRSS) introduced 
in our previous publications 57- 5^|: 



fa\ n 1 

U(r)=e( y -j +-e(l-tanh(fc (r- ( 7 1 ))), (1) 

where n = 14, fco = 10. a is "hard" -core diameter and 
<7\ = 1.35; 1.45; 1.55; 1.8 is soft-core diameter. 

In the remainder of this paper we use the dimensionless 
quantities: f = r/d, P = Pd 3 /e, V = V/Nd 3 = 1/p. As 
we will only use these reduced variables, we omit the 
tildes. 

In Refs. |59} it was shown that this system demon- 
strates anomalous behavior. A relation between phase di- 
agram and anomalous regions was also discussed in these 
articles. Our later publications gave detailed study of 
diffusion, density and structural anomalies in this sys- 



tem mO- 62 



We simulated the systems with four different step sizes 
o 1 and monitored the change in phase diagram and 
anomalous regions with increasing the step width. 

Importantly, in low temperature region the systems 
can demonstrate slow dynamics. In order to get reliable 



2 



results in the whole temperature range we used parallel 
tempering technic 64|. The system with o\ = 1.35 and 
a i = 1.45 were simulated at densities from p = 0.3 till 
p = 0.8 with step 8p = 0.05 (Sp = 0.025 in the vicinity of 
anomalous regions) and at temperatures from T = 0.12 
till T = 0.8. In the case of u\ = 1.55 the range of den- 
sities was p = 0.2 — 0.75. Finally, the density range of 
p = 0.1 - 0.8 was simulated for a\ — 1.8. Along every iso- 
chor we ran several parallel tempering runs with different 
temperatures. Lowe- Andersen thermostat was used dur- 
ing the equilibration time |65| ]. After equilibration the 
thermostat was switched off and the system evolved in 
NVE ensemble for production. After the production a 
trial change of temperature was made. Many data points 
were collected along every isochor (more then 100). Using 
these data points we constructed 9— th order polynomial 
approximation of internal energy, pressure and diffusion 
coefficient along isochors. These approximations were 
used in the following analysis. 

In order to compute the excess entropy of the liquids we 
used thermodynamic integration method 64|. We com- 
puted free energy along a high temperature isotherm by 
integrating equation of state of the system. After that we 
computed the free energies along isochors by integrating 
the U/T 2 function, where U is the internal energy 
Excess entropy was obtained as S e 



U-F el 
Nk B T 1 



64]. 



III. RESULTS 

Here we present the anomalous regions for the four 
systems studied. 



1.35 



The phase diagram and anomalies in this system were 
already reported in our previous publications 
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Here we summarize the previous results and place the 
anomalous regions on the phase diagram. We also think 
that it is necessary to repeat the results for this particular 
system for the sake of completeness. 

Figs. [1] (a)-(c) show the diffusion coefficient, pressure 
and excess entropy for the system with <j\ = 1.35. One 
can see that all three anomalies take place in the sys- 
tem. It is also evident that structural anomaly is more 
stable than the diffusion one since it disappears at higher 
temperatures. 

Fig. [2] places the anomalies regions on the phase dia- 



gram of the system reported in Refs. 57, 58]. One can 
see that all three anomalous regions appear exactly af- 
ter the low-density crystal bump. This allows to relate 
the anomalies to the phase diagram, i.e. the reason for 
the anomalies to appear can be the presence of the crys- 
tal bump inside the liquid phase. The principle question 
here is that liquids to the left from the bump and to the 



0.20 



0.15 



0.10 



0.05- 



0.00 



0.06 
0.04 
0.02 



v--- , 

. v 



J.3 0.4 0.5 0.6 0.7 0.8 



-■- 1=0.20 
-•- 1=0.25 

- 1=0.30 
-T- 1=0.35 

- 1=0.40 
-<- 1=0.50 



— i — i — i — 1 — i — 1 — i — 

0.3 0.4 0.5 0.6 



0.7 0.8 



4.0- 



3.5 



3.0- 



2.5 



2.0 




0.1 



0.2 



0.3 



0.4 



0.5 



0.6 




FIG. 1: (Color online), (a) Diffusion coefficient along a set 
of isotherms; (b) pressure along a set of isochors; (c) excess 
entropy along a set of isotherms for the system with ai — 1.35. 
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FIG. 2: (Color online). Location of anomalous regions at the 
phase diagram of the system with cri = 1.35. 



right have different local structures: while the highest ra- 
dial distribution function peak at the "left" liquid is the 
second one, corresponding to the soft core distance, at 
the "right" liquid it is the first peak located on the hard 
core distance. In between of these two liquids a region ap- 
pears where strong competition of two distances appears. 
It is this competition which brings to the appearance of 
anomalous behavior. 

The second point to notice is that the anomalous re- 
gions correspond to the picture proposed for water 
i.e. the diffusion anomaly region is inside the structural 
anomaly and the density anomaly is mainly inside the 
diffusion anomaly. 
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Figs. [3] (a)-(c) show the anomalies for the system with 
(Ji = 1.45. One can see from these pictures that the 
system behaves qualitatively similar to the previous case. 
However, the diffusion anomaly looks suppressed. At the 
temperatures as low as 0.15 the diffusion is close to a 
bend, but it is still monotonous. At the same time the 
density anomaly is rather pronounced in the system. 

The location of anomalous regions and the low density 
Face Centered Cubic (FCC) phase for o\ = 1.45 system is 
shown in Fig. [4] One can see that the diffusion anomaly 
has almost gone under the melting line. Only a small 
part of diffusion anomaly region is located in stable liquid 
phase. At the same time density anomaly is still very 
pronounced. It occupies a large region in the liquid part 
of the phase diagram. 
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FIG. 3: (Color online), (a) Diffusion coefficient along a set 
of isotherms; (b) pressure along a set of isochors; (c) excess 
entropy along a set of isotherms for the system with <ri = 1.45. 
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FIG. 4: (Color online). Location of anomalous regions at the 
phase diagram of the system with cri = 1.45. 



The most important fact extracted from Fig. @] is that 
the diffusion and density anomalies inverted with respect 
to each other, i.e. now diffusion anomaly region is inside 
the density anomaly one. The similar inversion of the 
order of anomalies is found in silica, however, in silica 
one can see the inversion of the order of structural and 
diffusion anomaly [HI . We can conclude that the case of 
unusual sequence of anomalies in silica is not unique, and 
the SRSS system is one more example of such unusual 
behavior. 



en = 1.55 



1.2- 
1.1- 
1.0- 
tj, 0.9- 
0.8- 
0.7- 
0.6 -. 



0.10 



(b) 



0.15 



0.20 



p=0.25 

p=0.275 

p=0.3 
p=0.325 



0.25 



In the case of u\ = 1.55 we do not find diffusion 
anomaly (Fig. [Sfa)). At low temperatures the diffusion 
along an isotherm develops a bend, but not a loop. At 
the same time both density and structural anomalies are 
present in the system (Figs.[5ljb)-(c)). Comparing it with 
the results for the previous two cases we conclude that 
the diffusion anomaly is completely under the melting 
line for the present case which makes it unobservable. 

From Fig. [5] we see that the density anomaly also ap- 
proaches the melting line. This allows to suggest that 
density anomaly also shrinks with increasing the step 
size. 
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In the case of o\ — 1.8 both diffusion and density 
anomalies are absent. Only excess entropy demonstrates 
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FIG. 5: (Color online), (a) Diffusion coefficient along a set 
of isotherms; (b) pressure along a set of isochors; (c) excess 
entropy along a set of isotherms for the system with o\ — 1.55. 
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FIG. 6: (Color online). Location of anomalous regions at the 
phase diagram of the system with ai = 1.55. 
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FIG. 8: (Color online). Location of structural anomaly region 
at the phase diagram of the system with a\ = 1.8. Note that 
the diffusion and density anomalies are absent for this case. 
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FIG. 7: (Color online). Excess entropy along a set of 
isotherms for the system with <ti = 1.8. 



anomalous increase along low-temperature isotherms 
(Fig. [7]). The situation with density anomaly looks ex- 
actly the same as with the diffusion one: the anomalous 
region shrinks with increasing the step and finally hides 
under the melting line. At the same time structural 
anomaly looks very stable with respect to the melting 
line: it is still located at the temperatures almost twice 
higher then the melting line maximum (Fig. [5]). 



IV. DISCUSSION 

Figs. O (a) and (b) summarize the results for all four 
systems. One can see that with increasing the step size 
(Ti the low density FCC phase and all anomalous re- 
gions move to lower densities and lower temperatures. 
At the same time the maximum temperature of both 
FCC phase and anomalous regions dramatically drops 
from tj i = 1.35 to o\ — 1.45 and decreases a bit from 
1.45 to 1.55. However, the rate of decreasing of different 
curves is different. The diffusion anomaly demonstrates 
the fastest decay with increasing the repulsive step size. 
As a result it quickly disappears under the melting line. 
The density anomaly is the second fastest to decay. As 
result it disappears somewhere between o\ = 1.55 and 
ai = 1.8. The structural anomaly is much more sta- 
ble than the diffusion and density ones and it does not 
demonstrate any precursors of disappearance at the mo- 
ment. It is interesting to note, that similar "resisting" 
behavior of the structural anomaly was found in Ref. [23j , 
where it was shown that making the soft-core of the po- 



tential from Refs. |25j, [26( steeper and more penetrable, 
the regions of density and diffusion anomalies contract in 
the T — p plane, while the region of structural anomaly is 
weakly affected. The authors of Ref. [23| concluded that 
a liquid can have anomalous structural behavior without 



having density or diffusion anomalies. However, in [23 1 
it was not found inversion of the order of anomalies with 
changing the parameters of the potential. It should be 
also noted that these results are consistent with results 
of Ref. [67j where it was proposed that for some of the 
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FIG. 9: (Color online). Location of anomalous regions at the 
phase diagram for (a) systems with <7i = 1.35 and o\ = 1.45; 
(b) o\ = 1.55 and ai — 1.8. 



core-softened systems anomalies go far under the melting 
line and become unobservable. 

In our previous work we proposed to consider the 
anomalous regions in scaled coordinates - p ■ af and 
T/T m ax, where T max is the temperature of maximum 
of the melting line of the low density FCC phase. This 
representation is very useful since it allows to see the rel- 
ative location of different regions in the phase diagram. 
Figs. [10] (a) and (b) give such representation for the sys- 
tems considered in the present paper. This diagram con- 
firms the conclusions of the previous paragraph. Interest- 
ingly, the relative height of the structural anomaly even 



increased from <j\ = 1.35 till 1.55. However, further in- 
crease to a\ = 1.8 pushes the structural anomaly down to 
the melting line. One can relate this qualitative change 
to the fact that the largest step involves in the interac- 
tion the second nearest neighbors which alters the system 
properties. Basing on this change we can not generalize 
the results to the step of arbitrary width: further increase 
of the step width will involve more and more neighbors 
into interaction which can induce more complex effects. 
In our previous publication we investigated liquid-liquid 
transition in core-softened systems in frames of pertur- 
bation theory (6|| . It was shown that the behavior of the 
system is "periodic" with increasing the step size: liquid- 
liquid transition takes place at some step size, then dis- 
appears and on further increase appears again. We can 
expect that the diffusion and density anomalies will oc- 
cur again at higher steps, however, this question requires 
further investigation. 



V. CONCLUSIONS 

In conclusion, this publication represents a detailed 
study of phase diagrams and anomalous behavior of 
Smooth Repulsive Shoulder System (SRSS) which be- 
longs to the class of core-softened systems. We inves- 
tigate the relation of the anomalous behavior with the 
phase diagram and the interaction potential parameters 
of the system. We show that both low density crystal 
phase and all anomalous regions shrink with increasing 
the step size. However, the rate of the decrease is dif- 
ferent. It is also shown that the order of the region of 
anomalous diffusion and the region of density anomaly is 
inverted with increasing the width of the repulsive shoul- 
der. 
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